Perturbations in stochastic inflation 



Kerstin E. Kunze ^ 

Departamento de Fhica Fundamental, 
and 

Institute Universitario de Fhica Fundamental y Matematicas (lUFFyM), 
I Universidad de Salamanca, 

^ ■ Plaza de la Merced s/n, E-37008 Salamanca, Spain 



o 



OO 

^ ■ Abstract 

^ i The perturbative approach to stochastic inflation is used to determine the spectrum of density 

fluctuations and gravitational waves due to the coarse grained field. The amplitude of the curvature 



in 



lO ' fluctuation spectrum, the spectral index and the running of the spectral index are in general found to 

. be smaller than in the standard approach to inflation. Furthermore, the amount of non-gaussianity 

' due to the second order perturbation is estimated. 

o ■ 
o 

' Inflation provides a mechanism in order to seed large scale structure formation. Quantum fluctua- 

■ tions of the inflaton are stretched beyond the horizon and are converted into classical perturbations. 

In the standard approach to calculating the spectrum of the resulting fluctuations expectation val- 
ues of the quantum fluctuations are identified with statistical averages of classical perturbations. 



1 Introduction 



;_j ■ Thus the amplitude of fluctuations is determined by the vacuum expectation value of the quantum 



field. The details of the quantum to classical transition are neglected. It is assumed that perturba- 
tions are classical instantaneously at horizon exit. The assumption of a classical behaviour outside 
the horizon has been found to be justified, as it is related to the particular behaviour of the mode 
functions in a de Sitter space-time (see for example PP). 

In the stochastic approach to inflation the dynamics of the quantum to classical transition is 
effectively described by a classical noise. The scalar field is coarse grained. Thus the field modes 
are split into a subhorizon part and a superhorizon part. The superhorizon part constitutes the 
classical, homogeneous coarse grained field driven by the stochastic noise due to the subhorizon 
modes leaving the horizon. As the inflationary universe keeps expanding rapidly more and more 
short wavelength modes are stretched beyond the horizon and subsequently contribute to the coarse 
grained, classical field. The coarse graining scale is at least of the order of the horizon. As shown 
in [2] in the case of a scalar field in a de Sitter space-time, the effective dynamics of the classical, 
coarse grained field is described by a Langevin equation with a white noise term. However, it 
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has been shown that the appearance of the white noise term is a direct consequence of the choice 
of window function to separate the sub- and superhorizon modes 0. Furthermore, models of 
stochastic inflation attempting to model in more detail the decoherence of the quantum field can 
be found, for example, in ||4j. 

Here, however, only the simplest model of stochastic inflation will be considered j2j. Thus the 
dynamics of the coarse grained field is determined by a white noise term. Recently, the equations 
of stochastic inflation for the coarse grained inflaton have been solved explicitly in a perturbative 
way (for earlier work in this context, see (6j). This allows to take into account as back reaction 
the evolution of the inflaton, namely the time variation of the Hubble parameter. In this work 
the perturbative approach is used to determine the spectrum of the resulting classical fluctuations 
due to the coarse grained field. Moreover, it is possible to estimate the amount of non-gaussianity, 
encoded in the parameter /tvl [Z|-^3i by including terms up to second order in the perturbative 
expansion. 



2 Curvature perturbation 

In the slow roll regime of stochastic inflation the evolution of the coarse grained scalar field (/? is 
described by the Langevin equation |2j, 



1 dV hI ^, , 



(2.1) 



where ^(i) is a white noise process with (^(t)) = and {(,{t)^{t')) = 6{t — t'). Furthermore, the 
Hubble parameter is given by H'^iif) = ^V{ip). Following 5 equation (HHJ) can be solved by 
the perturbative expansion 



(2.2) 



where ipc is the deterministic background value of the coarse grained field. Using the expansion 
1)2. 2|) in equation 1)2. 1|) results in j3] 
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(2.3) 
(2.4) 
(2.5) 



where a prime denotes the derivative with respect to (/7c- The solutions for the first and second 
order perturbations are given by 5^ 
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6ip2{t) 
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The spectrum induced by the first order fluctuation 6(fi can be calculated using the following 
expression for the curvature fluctuation, e.g. 



(2.8) 



Besides, since the curvature perturbation is constant outside the horizon it can be assumed that 
its spectrum is just a function of wave number. Thus 



Vn— 
q 



Ink 



(2.9) 



where k is some maximal scale, which corresponds to a minimal length scale that can be "probed" . 
This scale is determined by the coarse graining scale. The spectrum is calculated at the coarse 
graining scale, k = aaH, where cr is a parameter. Thus with dink = Hdt the spectrum is given by 



1 d 



in'), 



Hdt 

see for example Using the solution for 6ipi results in 

2 / TT\ 2 
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where (/Jq = V'c(O). Clearly, the first term is the standard result of the curvature perturbation 
spectrum due to a quantum field in a de Sitter background. The correction term is due to the 
evolution of the Hubble parameter in time. As can be shown, assuming H to be constant in time 
results in the correction term vanishing. In the standard approach the change due to the evolution 
of the inflaton is neglected. The variation in H is being considered only when deriving the spectral 
index and derived quantities of that (see for example PP). Moreover, choosing ipQ to be the value 
of the coarse grained field at horizon crossing, it is found that for the choice of the coarse graining 
scale to be of the order of the horizon there is no correction of the spectral amplitude. This is 
consistent with the picture that the spectrum in the standard approach is calculated at horizon 
crossing. In case that the coarse graining scale and the horizon are not the same equation ()2.11() 
implies that the spectrum is reduced with respect to the standard result. Modes with wavelengths 
between the horizon and the coarse graining scale are not taken into account. Thus the amplitude 
of the spectrum in the stochastic approach is in general reduced with respect to the amplitude in 
the standard approach. 

In the following the potential is assumed to be of the form 



Mp 



The mean square fluctuation at first order is found to be 



4^2 
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(2.13) 



In comparison, the amplitude of quantum fluctuations of a massless free scalar field in de Sitter 
space-time is given by Sep = H/2'k. However, it can be shown that at a time t = the mean 
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square first order fluctuation is given by {6(f1{H^^)) ~ Tliis is consistent witli the general 
picture that the stochastic differential equation (|2.1|) describes a random walk process with a step 
size of the order of H/2it over a characteristic time Here, this holds at first order in the 

expansion of the perturbation. 

The power spectrum of the curvature fluctuations Vn for the potential (|2.12() is found to be 



327r / Vq 




n+2 



This can also be written as V-r. = 1 ~ j ~ ^ } ^ (, A?^ j ' ^^^'^^ ^ = TB^ ('^) ^'^^ 

the slow roll parameters. In this case the spectrum depends explicitly on the value of the coarse 
grained field. Besides, the requirement V-ji > imposes an upper bound on the fraction 

implying the range 1 < < (l + n) ^' Thus this gives an upper bound on the size of the coarse 

graining domain. However, in principle, there should not be any limit. The upper cut-off is due to 
the fact that this is a perturbative approach and here only the first order is considered. Moreover, 
considering times of the order of the Hubble time leads to a change in the coarse grained scalar 
field of the order 

which is much smaller than the maximal bound imposed by the amplitude of the spectrum. The 
stochastic equation ()2.1|) is derived for de Sitter space-time. Using the stochastic approach to 
inflation to find the effect on the spectrum due to the evolution of the Hubble parameter can be 
interpreted as a perturbation around the de Sitter case. This is done varying the size of the coarse 
graining domain, that is varying the ratio 

For comparison, the spectrum in the standard approach is given by Vn = — 1 . Expression 



p 



1)2. 14() reduces to this for ( — I =1. The ratio of the amplitude of the curvature fluctuations in the 



<fic. 

stochastic approach over the one in the standard approach is shown in flguren(a). For (^^j > 1) 
the amplitude in the stochastic approach is always lower than in the standard approach. 
Whereas the spectral index given by 

d In Vt? 

is found to be 



^ ^ (n + 2)(n + 4)-n(n-2)(|^)S ^ X 



the running of the spectral index is given by 
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The quantities (|2.14|1 to (|2.18|) are calculated at = which is the field value A^^y e-folds before 
the end of inflation, when the characteristic scale of WMAP is of the order of the size of the coarse 
graining domain, ^p^r is given by 



The spectral index in the standard approach = 1 — 6e + 2r], with r/ = the second slow 

roll parameter, is given by 

ns,stand = 1 - ^(n + 2) (^^^ . (2.20) 

The spectral index calculated in the stochastic approach is shown in figure Q (b) for a range of 
values for N^r. The exact value of N]y depends on the history of the universe. Typically N]y 
is assumed to be the number of e-folds before the end of inflation when the present Hubble scale 
left the horizon, which should be a good approximation here. In ^1] it has been argued that the 
maximum range is given by 14 < N\y < 75. Other authors assume N\y to be in the range between 
50 and 60, see for example, ^Sj, or between 46 and 60 JHl- To see how values change in the case 
at hand we have chosen to be between 45 and 65. It is found that, for a given value of Nw, 
the larger the value ipo with respect to ^pw the smaller the spectral index. 

Assuming negligible running of the spectral index and negligible contributions from tensor 
perturbations the spectral index is found to be Ug = 0.951 ± 0.016 from the three year WMAP 
data only, = 0.948 ± 0.015 from WMAP data and data from the 2dF Galaxy Redshift Survey 
(2df) and Us = 0.948t|]:[}}^ from WMAP data and data from the Sloan Digital Sky Survey (SDSS). 
In general the spectral index in the stochastic approach depends on and Nw- It is always 
smaller than the one obtained in the standard approach, which only depends on N]y. This is true 
even at = 1. Furthermore, it is found that only the spectral index for a potential with n = 2 is 
compatible with observations. Moreover, in the case at hand the observational constraints can be 
satisified ( 
Nw = 65. 



satisified only for = 55 and N\y = 65, requiring < 1.08 for N\y = 55 and < 1.12 for 



The running of the spectral index in the standard approach ,li, dns^stand/dln k = 16ery — 24e^ — 
2,^^;, where = ^^y^^jr- is a third slow roll parameter, is given by 

dus^stand _ n^(n + 2) / ipc \ ~'^ ('2 9n 

d\nk ~ 327r2 \MpJ ' ^ ' ' 

The running of the spectral index in the stochastic approach is shown in figure ^c). Thus as in 
the standard approach, the running of the spectral index calculated in the stochastic approach is 
negative. For a chosen value of Nw, its modulus is increasing with increasing values of In 
figure 121 we have plotted in the case of a quadratic potential the running of the spectral index 
dus/dlnk versus the spectral index Ug in the stochastic as well as the standard approach. As 
can be seen in this figure the corresponding values for the running of the spectral index are much 
below observational bounds, assuming negligible contributions from tensor perturbations, = 
-0.055to:o3? from WMAP three year data only, ^ = -0.048 ± 0.027 from WMAP and 2df data 
and ^ = -0.060 ± 0.028 from WMAP and SDSS data [HI. Thus the assumption of negligible 
running of the spectral index is well justified. Similar results are found for the contribution from 
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Figure 1: (a) The ratio of the amphtude of the curvature perturbation due to the stochastic over 
the amphtude due to the standard approach „ . (b) The spectral index in the stochastic 

''^ stand 

approach, (c) The running of the spectral index in the stochastic approach. Graphs are shown for 
Nw = 45, 55, 65. Furthermore, n denotes the power of the potential. 



tensor perturbations (see next section) as can be seen from figure|31(b), which shows the tensor-to- 
scalar amplitude ratio versus the spectral index in the case of a quadratic potential in the standard 
and in the stochastic approach. Thus in order to satisfy observational bounds has to be close 
to 1. Therefore deviations from de Sitter space-time are small. 

3 Gravitational wave spectrum 

In this section the tensor-to-scalar amplitude ratio r is calculated which is another observational 
parameter. The gravitational wave amplitudes hij are effectively determined by two massless scalar 
fields ip^, ipx, respectively pp. In particular, hij = h^efj + hxe^j, where e"*" and are polarization 

tensors and /i+^x = ( 1 V'+,x • For simplicity, the index of V^+.x is omitted and the scalar fields 
are simply denoted by ^. Assuming that V is in the slow roll regime, at first order the following 
equation holds 



2tt 



-m, 



where the noise ^{t) is defined as before. This is solved by 



implying 



{^P\t)) = ^ 1^ drH^r). 
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Figure 2: The graph shows the running of the spectral index against the spectral index in the 
stochastic approach (lines) and the standard approach (points) for a quadratic potential for Nw = 
45,55,65. 



Thus the spectrum of gravitational waves, Vgrav = _ "^"'j^" ' , is found to be 

2 



16 / H \ 
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which equals the spectrum in the standard approach (cf. PP, e.g.). Hence the tensor-to-scalar 



amplitude ratio r = is found to be 



stock 



16e 



(3.26) 



Therefore, for — 1 the tensor-to-scalar amplitude ratio in the standard approach, Vgiand — 

(see for example ^01) is recovered. For all other values, r stock > r stand- This is due to the fact 
that in the stochastic approach the curvature spectrum amplitude is lower than in the standard 
approach. 

The tensor-to-scalar amplitude ratio calculated in the stochastic approach is shown in figure 
131(a) for different values of N\y- Assuming negligible running of the spectral index observational 
bounds on r are found to be: r < 0.65 from the three year WMAP data only, r < 0.38 from 
WMAP and 2df data and r < 0.30 from WMAP and SDSS data (17 (for the last value see also 
|16j). The constraint from the WMAP three year data only admits all three types of potentials 
considered in figure |31 (a) , up to a certain value of . The constraint from WMAP and 2df data 
marginally allows a quartic potential with N]^ = 45, which is already excluded by the WMAP 
plus SDSS data, independent of This stronger observational constraint can be satisfied up to 
a certain value of for potentials with n = A and Ny/ > 45 and for quadratic potentials. In 
figure 13 (b) the tensor-to-scalar amplitude ratio versus the spectral index is shown for a quadratic 
potential calculated in the stochastic as well as the standard approach. Values for r allowed by the 
observational constraints on the spectral index are of the same order in the stochastic and in the 
standard approach. 
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4 Non-gaussianity 



By construction, the first order perturbation 5ipi is gaussian with zero mean, {5'^i) = 0. Thus the 
curvature perturbation TZ is gaussian at first order. Non-gaussianity enters in TZ at second order in 
the expansion of 5ip (cf. equation H2.2() ). The nonUnearity can be estimated using the parameter 
fML defined by [Tj 



n = n„ 



ifNL (7^; 



(4.27) 



where TZg is a gaussian distribution with {IZg) = 0. Thus by construction (7^) = 0. In order to 
determine the non-gaussian contribution due to the second order perturbation, it is convenient 
to consider the curvature perturbation TZ = —M-5ip^ where dtp = Stpi + 5(p2 — {^f2)- Hence by 
construction (7^) = 0. Furthermore, the gaussian contribution is given by 



H 



(4.28) 



Thus ff^i^ is found to be 



NL 
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Taking into account only the dominant contribution in 5'p2-, coming from the second term in equa- 
tion ()2.7|) . yields to 



fNL 



This expression is exact for potentials with index n = 2 and n = 4, since in these cases H'" = 0. 
For a potential of the form V = Vq (^tup) equation 1)4. 3U() leads to, at Lpc = ^w, 
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The parameter | /atl | is shown in figure 01 (c) for several values of Nw and the potential index 
n. An estimate for (^-^^ is found using the expression for the power spectrum amplitude Vn (cf. 

equation 1)2. 14(1 ) together with Vtz = 2.4 x 10~^ As can be seen from figure|ni(c), /nl satisfies 
the bound from observations, —54 < /nl < 114 [T8] . 

For potentials with an 



For n = 4 the parameter I/nlI only depends very weakly on 



index n > 4 the parameter I/tvlI is increasing with increasing values of For n = 2 \fNL\ is 



a descreasing function for increasing values of 

^ > 1. 



The non-gaussian contribution requires that 
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Figure 3: (a) The tensor-to-scalar amplitude ratio r^toc/i in the stochastic approach, (b) The tensor- 
to-scalar amplitude ratio against the spectral index in the stochastic (lines) and in the standard 
approach (points) in the case of a quadratic potential, (c) The non-gaussianity parameter I/atl] in 
the stochastic approach. Curves are shown for Ny/ = 45, 55, 65. Furthermore, n denotes the power 
of the potential. 



These results can be compared with values for f^i found in the standard approach. Pertur- 
bation analysis up to second order in the standard approach to slow roll inflation leads to the 
estimate that f^L is of the order of the slow roll parameters |Hl-d|- Inl is found to be 
/iVL ~ n^""* where ret is the spectral index of the gravitational wave spectrum and 

< / < I characterizes the configuration of the momenta which enters into the calculation of 
the three-point correlation function of the curvature perturbation. Choosing /tvl = | (equilat- 
eral triangle) results for the range of parameters under consideration here {Nw = 45, 55, 65 and 
re = 2, 4, 6) in a non-linearity parameter /tvl ~ C(10~^). Thus it is of the same order of magnitude 
as fNL calculated in the stochastic approach (cf. figure El (c)). 

5 Conclusions 

Recently, the equations of stochastic infiation have been solved in a perturbative way [3]. This 
allows to take into account the effect on the perturbations due to the evolving Hubble parameter. 
Here, the perturbative approach to stochastic inflation has been used to determine the spectrum 
of the curvature perturbation due to the coarse grained inflaton as well as the gravitational wave 
spectrum. 

The spectrum of curvature fiuctations has been calculated using the solution for the perturbation 
of the coarse grained field up to first order in the noise. The resulting expression for the spectrum 
contains an additional parameter in the form of a value 999 of the coarse grained scalar field (/?. This 
has been interpreted as the value of the deterministic coarse grained inflaton tpc at the horizon. 
This interpretation is motivated by the fact that the resulting perturbation spectrum reduces to 
the standard expression in the case = V^o- In the standard approach the curvature perturbation 
spectrum is found using the spectrum of vacuum fiuctations in a de Sitter background calculated at 
the horizon. The coarse graining of the scalar field in the stochastic approach to inflation introduces 
a coarse graining scale. This is at least of the order of the horizon. The spectrum is calculated 
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at the coarse graining scale. In general the amplitude of the spectrum of curvature fluctuations 

calculated using the stochastic approach is reduced, in comparison with the standard approach, 
for ^ > 1. This results from the horizon and the coarse graining scale not being the same in this 
case. Modes with wavelengths between the horizon and the coarse graining scale do not contribute 
to the spectrum. 

The expression for the spectral index depends in general on the parameter ^pQ. Even at (pc = (fo 
the spectral index is smaller than the one obtained in the standard approach. Moreover, only for 
polynomial potentials with a potential index n = 2 observational bounds can be satisfied. The 
running of the spectral index calculated in the stochastic approach is negative as it is the case in 
the standard approach. Apart from its dependence on the index of the potential and the number 
of e-folds before the end of inflation, at which the coarse grained field is evaluated, its value also 
depends on 

The tensor-to-scalar amplitude ratio is in general larger in the stochastic than in the standard 
approach for ^ > 1. This is due to the fact that whereas the gravitational wave amplitude is 
found to be the same as in the standard approach the scalar perturbation amplitude is lower in the 
stochastic approach. 

Finally, the non-gaussian contribution due to the second order perturbation in the stochastic 
approach to inflation has been calculated. By construction perturbations at first order in the noise 
are gaussian since the noise is gaussian. Non-gaussianity enters at second order in the noise. The 
resulting expression for the parameter {/nlI reqires ^ > 1. For potentials with a potential index 
n = 2 it is a decreasing function with increasing for n = 4 it is only very weakly depending on 
the value of ^ and for n > 4 it is an increasing function with The resulting values for I/nl] 
are compatible with bounds from observations. 
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